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Cerium 4/ electronic spin dynamics in single crystals of the heavy-fermion system CeFePO is 
studied by means of ac-susceptibility, specific heat and muon-spin relaxation (/iSR). Short-range 
static magnetism occurs below the freezing temperature T g w 0.7 K, which prevents the system 
from accessing the putative ferromagnetic quantum critical point. In the /iSR, the sample-averaged 
muon asymmetry function is dominated by strongly inhomogeneous spin fluctuations below 10 K 
and exhibits a characteristic time-field scaling relation expected from glassy spin dynamics, strongly 
evidencing cooperative and critical spin fluctuations. The overall behavior can be ascribed neither 
to canonical spin glasses nor other disorder-driven mechanisms. 

PACS numbers: 71.27.+a, 64.70. Tg, 76.75.-K, 75.50.Lk 



A long-standing question in the field of quantum criti- 
cality is whether a ferromagnetic (FM) quantum critical 
point (QCP) generally exists and, if not, which are the 
possible ground states of matter that replace it. Quan- 
tum critical points occur when a material is continuously 
tuned with an external parameter (pressure, magnetic 
field, etc.) between competing ground states at zero tem- 
perature [H, 0). An FM-QCP then exists when it is pos- 
sible to shift the Curie transition temperature Tc of a 
ferromagnet continuously to zero where a second order 
quantum phase transition takes place. Quantum phase 
transitions occur at zero entropy and are driven by quan- 
tum rather than thermal fluctuations. These fluctuations 
diverge at the QCP modifying the excitation spectrum of 
a metal and leading to a fundamental instability of Lan- 
dau's Fermi liquid (FL) Q. Typical signatures of such 
a behavior are observed in magnetic, thermal and trans- 
port properties and are referred to as non-Fermi-liquid 
(NFL) phenomena 

Although there is clear evidence for the existence of 
antiferromagnetic (AFM) QCPs, the FM-QCP case is 
controversial. In recent years, substantial experimen- 
tal and theoretical efforts were made to further inves- 
tigate this problem. However, a wide range of possi- 
bilities exists. On theoretical grounds, a 3-dimensional 
(3D) FM-QCP is believed to be inherently unstable, ei- 
ther towards a first order phase transition or towards 
an inhomogeneous magnetic phase (modulated/textured 
structures) l5H7| . Similar results have been obtained in 
2D 0, 0, Several clean (stoichiometric) magnetic 
transition-metal compounds, like MnSi OH or ZrZn 2 [IH, 



show NFL behavior close to a FM instability, but the 
transition changes into a first order one. In other sys- 
tems the existence of a FM-QCP has been proposed, 
most notably in Nbi_ y Fe2+ y |12|, Zri-^Nb^Z^ [l3( or 
SrCo^Gei-^P^ [3 where the FM-QCP is attained by 



chemical substitution. However, in these cases the influ- 
ence of disorder remains ambiguous. 

More appropriate candidates for the study of FM- 
QCPs are U-, Yb or Ce-based /-electron metals @, 0], 
since in these materials the NFL signatures are much 
more pronounced due to their heavy-fermion (HF) char- 
acter. However, while there is quite a number of U- 
based systems showing either a first order FM transi- 
tion (UGc 2 , UCoAl, UCoGe) or indications for a FM 
QCP (UCug^Pdz 0, URhi-ajRusGe, 0), the num- 
ber of Yb-based systems close to a FM QCP is very lim- 
ited (YbNi 4 P 2 [ill, YbCu 2 Si2_p|). Several systems, like 
CeRu 2 Ge 2 or CeRuPO dJwhere the FM transition 
temperature is suppressed to T = by hydrostatic pres- 
sure, exhibit a change into AFM order before reaching 
the QCP. There are Ce-based alloys (CePdi-aRh^ (22J) 
and also d-electron metals (Nii-^Vj; [23[) where it seems 
that local disorder-driven mechanisms such as Kondo dis- 
order or the quantum Griffiths phase (QGP) scenario are 
responsible for the NFL properties [24 - 1261 ] . Broad and 
strongly T dependent NMR and ^SR linewidths are in- 
dicative for such disorder-driven mechanisms. As a con- 
sequence spin-glass-like behavior is often found, e.g., in 
CcPdi-^Rhj;, and power-law corrections to the thermo- 
dynamic and transport properties as well as in the local 
spin dynamics are expected in a broad region across the 
putative QCP. The global phase transition then becomes 
smeared [ll, [27|, Sil • 

In this context, the layered Kondo-lattice system Ce- 
FePO is a unique candidate for studying FM-QCPs, since 
it is a clean non-magnetic (non-supcrconducting) HF 
metal located very close to a FM-QCP with strong FM 
fluctuations [2^, [30]. CeFePO is a homologue of the qua- 
ternary iron pnictides. It evolves from a long-range or- 
dered FM ground state, when a small amount of arsenic 
is substituted for phosphorus 31-33J. Less As concen- 
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Figure 1: (Color online) (a) Ac-susceptibility x'(^) OI the PC 
sample with its imaginary part x" {T) (inset), (b) x'(T) vs 
T for the SC sample with B 1 c. (c) x'(T) vs T for the SC 
sample with B\\c. (d) Specific heat of the PC sample plotted 
as C/T vs T. 



tration leads to a continuous decrease of Tq culminat- 
ing into a putative FM-QCP. The isovalent As substi- 
tution not only introduces a volume effect (shortening 
mostly the c-axis) but also increases locally the hybridiza- 
tion strength between the trivalent Ce-4/ and the Fe-3(i 
conduction electrons leading to an enhancement of the 
Kondo temperature, which is approximately 1 K for Ce- 
FePO (HI. Below this temperature, the susceptibility 
and the Knight shift become field dependent and the 
NMR line width broadens. This identifies the onset of 
short-range FM correlations essentially within the basal 
plane, evidencing a strong anisotropy (2^ |. which could 
also be confirmed by recent NMR measurements on ori- 
ented powder [35$ . The ground state of CeFcPO was 
found to be a paramagnetic heavy FL (29j . 

In this letter, we present a comprehensive study of 
the ac-susceptibility (%'), specific heat (C) and muon- 
spin relaxation (/xSR) measurements on recently grown 
high-quality single crystals of CeFePO. We find evidence 
of strongly inhomogeneous spin fluctuations starting be- 
low 10 K and of a spin-glass-likc freezing at T g w 0.7 K 
which prevents the system from accessing the putative 
FM-QCP. The observed time-field scaling of the muon 
asymmetry points to a cooperative mechanism and to 
the presence of critical spin fluctuations. The overall be- 
havior can not be ascribed to either canonical spin glasses 
or to other disorder-driven mechanisms. The physics of 
CeFePO is different from other candidate systems where 




Figure 2: (Color online) Frequency dependence of x'{T) for 
the PC sample. Inset: Arrhenius plot for the PC sample 
(black) and for the SC sample (green) with B _L c. 



the putative FM-QCP is avoided by a first order phase 
transition or a transition into an AFM ordered state. 

The samples were synthesized by means of a two-step 
Sn-flux method. The small crystals were powdered and 
pressed into pellets which are referred to as polycrys- 
talline (PC) sample in the subsequent discussion. The 
large crystals were oriented and glued together with silver 
paint to form a larger "single crystal" (SC). X-ray pow- 
der diffraction confirms the ZrCuSiAs-structure type [3l[ . 
The synthesis conditions were different from the samples 
investigated previously pil 31 1. However, 



the samples 

PC and SC as well as the sample of Ref. [29[ are in- 
distinguishable within the resolution of energy dispersive 
X-ray (EDX) measurements which confirm the stoichio- 
metric ratio 1:1:1:1. Low temperature x'(T,B), C(T,B) 
and //SR(T, B) were measured in 3 He- 4 He-dilution refrig- 
erators. A commercial SQUID-VSM (Quantum Design) 
was used to measure x' above 1.8 K. The /xSR experi- 
ments were performed on the 7rM3 beam line at the Swiss 
Muon Source at the Paul-Scherrer-Institut, Switzerland. 
The fiSK measurements were executed in zero magnetic 
field and in applied magnetic fields up to 0.75 T parallel 
to the initial muon-spin polarization (LF-/iSR). 

The first evidence of spin freezing is seen in the T de- 
pendence of x'(T) for the PC sample (Fig.QJi). At B = 
a distinct peak is found at T g = 0.67(1) K where the sus- 
ceptibility reaches values as high as 2.8 • 10~ 6 m 3 /mol. 
The x' (T) value at 2 K is three times larger than the one 
measured in the sample of Ref. [2^] . With increasing field 
its amplitude decreases and T g shifts slightly and, above 
0.8 T, x'(T) flattens. The small hump between 3 and 
9K is only seen at B = in the PC sample. It is most 
likely due to an impurity phase, as it is completely sup- 
pressed in a small field of 0.02 T while T g does not change. 
Moreover, a comparable feature is absent in the SC case. 
Dissipativc effects are corroborated by a peak in x"(^) 
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at about 0.45 K (inset of Fig. [T^,), which is a significantly 
lower temperature than that of the maximum in x'(T). 
The same effects and B dependence are observed in the 
SC sample, but at a higher T g = 0.92(2) K. This could be 
due to a very tiny difference in stoichiometry, since both 
the PC and SC samples were taken from small and large 
crystals, respectively, of the same batch which likely form 
at different times during the growth. The susceptibility 
is very anisotropic (x± c /x\\c ~ 5), and with B _L c it 
reaches a high peak value of 9 • 10 -6 m 3 /mol (Fig. [TJd and 
c) which is much larger than in the PC case. This ob- 
servation confirms the presence of anisotropic FM spin 
fluctuations, which are much stronger along the basal 
planes. To check the bulk nature of the freezing, we have 
measured the specific heat of the PC sample (Fig. QJ1). 
A broad maximum emerges at about 0.55 K in a C/T 
vs T plot. While at 0.1 T the maximum is unchanged, 
larger fields suppress it. The entropy difference between 
the zero-field and the high-field curves is small, about 1% 
of i?ln(2), which is due to the Kondo screening of the Ce 
moments. The data of the polycrystalline sample inves- 
tigated in Ref. [29| (gray data points in Fig. [TJi) resemble 
the C/T behavior at high temperatures, while below 3 K, 
C /T remains constant without any indication of freezing. 

We have performed frequency (/) dependent measure- 
ments at B = on the PC and SC samples to investigate 
the spin dynamics (Fig. [2]). Similar to spin glasses, the 
maximum in x'(^) shifts to higher temperatures as the 
excitation frequency is increased, while its amplitude de- 
creases. T g and / are shown in an Arrhenius plot in the 
inset of Fig. ^ from which the frequency shift 5 [36[ can 
be evaluated. For the PC and the SC sample we ob- 
tain Spc = 0.085(11) and dsc = 0.065(16), respectively. 
These values are larger than those found for canonical 
spin glasses (5 w 0.005. ..0.06), yet they are below typi- 
cal values of superparamagncts (S w 0.3) [36|. A qual- 
itatively similar behavior was found for CePdi_ x Rh x , 
where S increases when approaching the critical point 
at x = 0.87 [22j | . It is interesting to note that in 
CeFePO the sample with lower T g exhibits larger 8, 
too. This suggests that the QGP scenario might be 
also applicable in CeFePO. On the other hand, the mag- 
netic anisotropy and the lack of evidence for FM cluster 
formation in zero-field-cooled/ficld-cooled magnetization 
measurements (not shown) rule out such a mechanism 
in CeFePO. Moreover, we could not fit our data with 
X oc C/T oc T A_1 with < A < 1 in any reasonable 
T-range [IBj]. 

To study the microscopic nature of the low-T mag- 
netism in CeFePO, /iSR measurements in zero and small 
longitudinal fields were performed on the PC sample. 
/iSR in small fields is dominated by thermally excited 
Ce-4/ electronic spin fluctuations that couple to the im- 
planted muons. Figure [3ja, displays the temperature evo- 
lution of the normalized muon-spin asymmetry function 
G(t, B) at a constant field of B = 0.01 T which is suffi- 
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Figure 3: (Color online) (a) Normalized muon-spin asymme- 
try function G(t,B) at an applied field B — 0.01 T for repre- 
sentative temperatures above and below T g . Solid lines are fits 
according to Eq. [TJ For clarity, the different curves are shifted 
subsequently by 0.33. (b) T dependence of the dynamic ^tSR 
rate Xl- The blue arrow marks the magnetic transition from 
dynamic to static magnetism at T g — 0.70(3) K. (c) T de- 
pendence of the static /xSR rate At- (d) T dependence of 
the exponent j3 in Eq. [TJ The dotted line denotes the value 
/3 = 0.5. 



cient to quench the weak static relaxation due to nuclear 
dipole fields at the muon sites, leaving only the dynamic 
and static contributions due to the electronic Ce-4/ mag- 
netic moments. The quantitative analysis takes into ac- 
count both static and dynamical fields: The static relax- 
ation dominates G{t 1 B) at short times t, while at long t 
the relaxation rate probes only the dynamic spin fluctu- 
ations, i.e., the Fourier transform (FT) of the dynamic 
spin-spin autocorrelation function q(t) = (Si(t) ■ S,(0)). 
To account for both contributions we use the following 
fitting function 



G(t, B) = G 1 exp[-(A T *)] + G 2 cxp[-(A L £)' 3 ] 



(1) 



with static (transversal) and dynamic (longitudinal) re- 
laxation rates At and Al, respectively. The fits provide 
a very good description of the experimental data (solid 
lines in Fig. [3Ji). At low T the spectra show nearly no 
muon-spin relaxation at long times, i.e., very small re- 
laxation rates Xl- Upon increasing T, Xl increases and 
reaches a maximum at T g = 0.70(3) K (Fig. (3}d) in agree- 
ment with T g found in x'{T). Subsequently, it decays 
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Figure 4: (Color online) (a) Field dependence of the muon- 
spin asymmetry function G(t, B) in CeFePO at T = 0.02 K 
and (b) T = 0.9 K. (c) The same data as in (b), plotted as a 
function of the scaling variable t/B ' 5 . 



following a T _1 behavior up to 10 K after which no dy- 
namic relaxation is observed as expected from/in accor- 
dance NMR and susceptibility experiments [29| . On the 
contrary, the static component At increases steeply be- 
low T g up to a value of 7 /is -1 (Fig. [3p). Its behavior 
resembles that of the magnetic order parameter when 
entering a magnetically ordered phase. The /3 value of 
about 0.5 for T > T g (Fig. [3Ji) indicates a broad in- 
homogeneous distribution of fluctuating dynamical local 
fields (or relaxation rates) [13 • In the ordered phase (3 
increases, reaching a value of about 1.7 at T = 0.02 K, in- 
dicating that the spin fluctuations become static. Figure 
Ha displays G(t, B) at 0.02 K. At B = 0, the absence of 
a spontaneous muon-spin precession frequency indicates 
short-range magnetic order. Such strongly damped /iSR 
signal allows for an estimation of the magnetic coherence 
length £ < 10 • a, where a is the lattice constant [38j [ . 
The observation of a 2/3 and 1/3 signal fraction below T g 
proves that 100 % of the sample volume shows static mag- 
netic order (39j . The muon-spin relaxation is completely 
suppressed at B = 0.75 T demonstrating that the inter- 
nal field distribution is static in nature at T = 0.02 K. 
Increasing T, a dynamic contribution to the muon-spin 
relaxation develops (Fig. Hp) . 

Glassy spin dynamics generally result in long-time corre- 
lations with distinct signatures when T g is approached 
from high T [4(j. Theoretically, q(t) is predicted to 
exhibit power-law q(t) = ct~ a or stretched exponen- 
tial q(t) = cexp[-(At) K ] behavior at T > T g , that in 
both cases can lead to a characteristic time-field scal- 
ing G(t,B) = G(t/B 7 ) after Fourier transforming q(t), 
where 7 < 1 and 7 > 1 for power-law and stretched ex- 
ponential correlations, respectively j4jj. If this equation 
is obeyed a plot of G(t, B) versus t/B"' at T > T g will be 
universal. FigurcHp displays G(t, B) at T = 0.9 K, which 



is slightly above T g , both in zero field and in magnetic 
longitudinal field between 0.01 and 0.75 T. The observed 
B dependence corresponds to a measurement of the FT of 
q(t) over the frequency range ^^B/2-k rj 1.4 — 100 MHz, 
where 7 M = 2-7T x 135.53 MHz/T is the muon gyromag- 
netic ratio. As shown, the relaxation slows with increas- 
ing B. For low enough B, the B dependence is expected 
to be due to the change of j^B rather than an effect of 
field on q(t) j3?| • A breakdown of time-field scaling is 
expected for high fields where q(t) is directly affected by 
the applied fields. Figure Up shows the same muon-spin 
asymmetry data as a function of t/B 1 . For 7 = 0.5(1) 
the data scale well over nearly 2.5 orders of magnitude 
in t/B 1 for all applied fields except for 0.75 T, as ex- 
pected for fields with [IbB > ksT, which should affect 
q(t). The obtained scaling exponent 7 = 0.5(1) < 1 
implies that, within the fjSK frequency range, q(t) is 
well approximated by a power law, suggesting cooper- 
ative and critical spin fluctuations rather than a distri- 
bution of local fluctuation rates (37J • This is in contrast 
to j3 = 0.5 which indicates a broad inhomogeneous dis- 
tribution of local fluctuation rates. The cooperative be- 
havior is supported by the fact that CeFePO is a stoi- 
chiometric system and the narrow NMR linewidth proves 
that it is locally not disordered (29|. Short-range cor- 
relations set in below 10 K, broadening the linewidth, 
in agreement with the //SR results. The value of 7 
seems to be weakly T dependent (it is 0.4 at about 2K, 
not shown) which suggests slow quantum rather than 
thermal fluctuations. From the B dependence of Al, 
the spin autocorrelation time t c can be estimated us- 
ing X L (B) = (2 7 ^ 2 > c )/[l + {llB 2 T 2 c Y] 0. Here, 
B\ oc (t) describes the time- varying local magnetic field 
at the muon site due to fluctuations of neighboring Ce- 
4f moments, with a local time averaged second moment 
7 M (_B[ 2 oc ). For fojj ^> ksT, the fluctuation-dissipation the- 
orem relates r c to the imaginary component of the lo- 
cal q- independent /-electron dynamic susceptibility, i.e., 
t c (B) w (k B T/^i 2 B ) (x"(w)/w) |4i|. A fit to the data (not 
shown) yields r c w 1.8 x 10 _8 s and p = 0.67. The value 
for t c indicates very slow or glass- type spin dynamics. 

In conclusion, we have shown that single crystals of Ce- 
FePO, located close to a FM instability, show spin-glass- 
like freezing. It is evidenced by the frequency-dependent 
peak in x'(^) & t T g , the broad maximum in C(T)/T 
as well as clear signatures in G(t,B), which clearly evi- 
dences the transition from dynamic to short-range static 
magnetism. The frequency shift of T g suggests values 
slightly larger than for canonical spin glasses, indicating 
the presence of large fluctuating regions above T g , and 
the time-field scaling strongly suggests cooperative be- 
havior. Our results imply that the putative FM-QCP is 
avoided in a new manner: We do not observe a first order 
phase transition or AFM order, but rather a transition 
into a short-range ordered state. Moreover, the magnetic 
anisotropy, the lack of evidence for FM clusters and the 
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time-field scaling rule out a disorder-driven scenario (e.g., 
the QGP) as the mechanism underlying the spin dynam- 
ics in CcFcPO. We might have in CcFcPO a possible 
combination of both scenarios: The close proximity of 
CeFePO to a FM instability and its magnetic anisotropy 
seem to drive the system to develop short-range magnetic 
correlations which might have their origin in the mecha- 
nism described in Refs. 0-0- Below T g , magnetic short- 
range order then forms with a certain texture (e.g., homo- 
geneous cobbled magnetically ordered regions of different 
sizes) which would explain the spatially distributed /xSR 
rates. 
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